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ANODE  PLASMA  EXPANSION  IN  PINCH-REFLEX  DIODES 


Introduction 

Anode  plasma  dynamics  is  important  in  all  pulsed-power  ion  diodes.  In 
most  diodes,  it  leads  to  anode-cathode  (A-K)  gap  closure  and  to  the  ion  pulse 
termination.  It  may  also  have  undesirable  effects  on  ion  beam  brightness 
since  ions  extracted  from  a  non-uniform  surface  will  not  be  accelerated  in  the 

A-K  gap  for  the  same  amount  of  time  and  not  necessarily  in  the  same 

1-4 

direction.  Few  studies  ,  experimental  or  theoretical,  have  been  made  of 
this  critical  phenomenon.  One  study^  has  determined  some  experimental 
characteristics  of  this  expansion.  First,  large  closure  velocities  up  to 
30  cm/ysec  have  been  observed.  Second,  this  large  expansion  is  observed  only 
late  into  the  pulse  (about  30  nsec  into  the  ion  pulse)  and  finally  it  takes 
place  preferentially  near  the  anode  axis  of  symmetry.  The  30  cm/usec 
measured  would  correspond  to  a  thermal  expansion  of  a  0.4  keV  CH2  plasma  which 
cannot  be  obtained  by  ohmic  heating  alone  for  the  measured  plasma 
parameters.  The  two  main  theoretical  attempts  at  explaining  the  anode  plasma 
expansion  have  not  been  very  conclusive.  Nardi  et  al^  have  shown  that 
10  cm/jsec  expansion  takes  place  near  the  anode  center  because  of  the  electron 
energy  deposition  in  that  region  but  they  did  not  reproduce  the  larger 
velocities  observed  nor  the  long  delay  in  the  expansion.  Prono  et  al  have 
studied  a  charge-exchange  mechanism  to  account  for  the  large  expansion 
velocity  but  this  mechanism  does  not  explain  the  favored  central  expansion  nor 
the  large  delay  in  the  expansion.  In  the  present  work,  we  describe  a  model 
based  on  the  dominant  role  played  by  the  (jxB)  force  in  the  anode  plasma 
expansion.  In  particular,  the  main  results  from  this  model  are  consistent 
with  large  final  expansion  velocities,  with  large  expansion  near  the  anode 
Mamucript  approved  June  20,  1983, 
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center  and  with  a  long  delay  before 
and  the  assunptions  are  presented. 


fast  expansion  sets  in.  First,  the  nodal 
Then  results  and  discussions  of  the  nost 


influential  parameters  follow 


Model 


The  geometry  of  the  model  is  shown  in  Fig.  1.  The  anode  foil  is  planar 
and  made  of  a  non-conducting  material  (polyethylene).  When  electrons  from  the 
cathode  accumulate  on  the  anode  surface,  flash-overs  occur  and  an  anode  plasma 
is  formed  just  in  front  of  the  anode  foil.  This  anode  plasma  serves  as  the 
source  of  the  ions.  When  ions  leave  the  plasma  generating  an  ion  current  in 
the  diode,  the  thermal  electrons  that  are  left  behind  flow  radially  inwards 
through  the  anode  plasma  where  they  are  collected  on  a  conducting  rod  at  r 
*  rmin'  The  rest  the  diode  current  flows  from  the  cathode  in  the  form  of  a 
relativistic  electron  beam  (REB)  which  is  assumed  to  pinch  in  vacuum  from  the 
diode  radius  R  to  a  radius  comparable  to  tm^n*  For  z  <  0,  all  the  diode 
current  is  carried  in  the  conducting  rod  and  therefore,  the  magnetic  field  due 
to  the  total  current  develops  behind  the  anode  plasma.  On  the  other  hand,  for 
z  >  0,  only  a  fraction  of  the  total  current  contributes  to  the  magnetic 
field.  Because  of  that  configuration,  a  gradient  in  magnetic  pressure 
develops  in  the  plasma.  The  force  on  the  plasma  due  to  this  gradient  is 
described  by  the  local  (j  x  B)  where  j  is  due  to  the  thermal  electron  current 
in  the  anode  plamsa.  Notice  that  in  the  present  model,  the  REB  does  not 
deposit  its  energy  in  the  anode  plasma  but  does  contribute  to  the  magnetic 
field  there.  The  anode  plasma  is  collision-dominated  and  thus  a  fluid 
description  is  used  to  study  its  motion  in  the  z-direction.  For  simplicity, 
we  have  performed  one-dimensional  studies  of  the  anode  plasma  motion  at 
various  radii,  approximating  a  two-dimensional  picture  by  assuming  that  the 
thermal  electron  current  has  only  a  radial  component. 

The  equations  which  describe  the  anode  plasma  are  the  magnetohydrodynamic 
equations^  which  express  conservation  of  mass,  momentum  and  energy.  A 
magnetic  diffusion  equation  closes  the  system.  In  the  momentum  aquation,  both 
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ANODE  PLASMA 


Fig.  1.  Schematic  of  the  anode  plasma  model  for  pinch-reflex  diode-plasma. 
Expansion  is  studied  in  the  z-direction  as  a  function  of  radius. 
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thermal  pressure  gradient  and  jxB  act  on  the  plasma.  In  the  energy  equation, 
radiation  losses  have  been  neglected  and  the  energy  input  into  the  plasma 
comes  from  the  Ohmic  heating  term  as  well  as  jxB  acceleration.  The  initial 
and  boundary  conditions  which  complete  the  set  of  assumptions  are  as 

follows.  A  low  temperature  is  assumed  for  the  C2H4  plasma  at  t  ■  0  (T  ■  0.5 

eV)  as  well  as  near-solid  density  (n^  *  nsolid^^  which  is  assumed  uniform 
over  a  distance  of  10  ^  cm  and  corresponds  roughly  to  an  energy  deposition^  of 
0.1  J/cm  .  The  electron  and  ion  current  waveforms  are  based  on  ref.  4.  They 
rise  linearly  in  time  in  30  nsec  and  remain  constant  thereafter  as  shown  in 
Fig.  2.  As  for  boundary  conditions,  the  ion  current  density  jz  for  z  >  zy 
where  zv  denotes  the  plasma-beam  interface  is  assumed  to  be  inversely 
proportional  to  r.  This  relationship  has  been  demonstrated  both 
experimentally  and  in  particle  simulation  .  The  resulting  B  field  at  z  »  zv 

due  to  the  ion  current  is  independent  of  radius.  The  B  field  due  to  the  REB 

current  is  inversely  proportional  to  r  for  z  -  0  and  z  *  Because  of  the 

time  dependence  of  the  currents  (Fig.  2),  the  B  field  increases  with  time  both 
at  z  *  0  and  z  ■  zv  and  is  larger  at  z  ■  0.  The  boundary  condition  to 
determine  the  position  of  zy  follows  from  both  plasma  expansion  and  plasma 
erosion.  The  expansion  is  determined  simply  by  the  plasma  boundary  velocity 
v  <zv).  The  anode  plasma  erosion  is  due  to  the  ion  beam  extraction  by  the 
diode  electric  field  at  z  >  zy  which  prevents  plasma  electrons  from  following 
the  ions.  The  expression  for  zy  is  then 


v 


F 


dz 

v 

dt 


v  (z  ) 

p  V 


£i 

ni 


where  the  subscript  i  refers  to  the  ions.  From  that  expression,  and  for 

2  7 

typical  peak  parameters;  -  100  kA/cm  ,  vp  =  10'  ca/sec,  the  plasma  number 
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density  at  zv  is  n^  »  10^  cm"^  or  p  s  10-^  g/cm^  for  a  CH2  plasma. 

Other  features  of  the  model  concern  its  treatment  of  the  equation  of 
state  and  of  the  electrical  resistivity.  The  equation  of  state  for  CH2  was 
taken  from  standard  SESAME  tables  except  that  at  low  temperatures  and  high 
density,  corrections  were  made  to  take  into  account  the  fact  that 
dissociated  into  CH2,  then  into  C  +  H  +  H.  Electrical  resistivity  was  assumed 
to  be  classical  and  was  supplemented  by  electron-neutral  collision  resistivity 
at  low  temperatures  and  by  anomalous  resistivity  at  high  temperature  and  low 
density. 
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Results 


Results  were  obtained  for  1^  «  I  <*  ’  >0  kA,  t  (linear  ramp)  *  30 

max  emax  P 

nsec  (see  Fig.  2),  R  a  5  cm  and  initial  plasma  mass  density/unit  area  *  3.5  x 

—6  2 

10  g/cst  .  Plasma  expansion  was  calculated  for  the  different  radii  r  =  1, 

2.5  and  4  cm.  The  following  is  a  brief  summary  and  analysis  of  the  results. 

In  all  cases,  pressure  gradients  dominated  the  j  x  B  force  in  the  first  few 
nanoseconds  and  gave  rise  to  expansion  velocities  of  a  few  times  the  thermal 
sound  speed  or  a  few  cm/ysec. 

When  the  j  x  S  force  becomes  important,  it  first  has  the  effect  of 

slowing  down  the  plasma  in  the  vicinity  of  zv*  In  order  to  explain  that 

effect,  we  use  as  an  example  the  information  in  Fig.  3  where  density, 

velocity,  8-field  and  current  density  profiles  are  shown  as  a  function  of  z  at 

t  *  14  nsec  and  r  *  2.5  cm.  The  main  feature  of  the  8-field  profile  in  Fig. 

c  3B 

3a  is  the  dip  inside  the  plasma.  Because  —  ,  we  see  that  for 

z  <  Zj,  (j  x  B)  accelerates  the  plasma  outwards  whereas  for  z  >  z^,  (j  x  3) 
slows  down  the  plasma.  The  reasons  for  this  dip  are  the  following:  near  z  * 
0,  the  plasma  density  Is  very  high,  the  temperature  rises  very  slowly  and  the 
magnetic  field  diffuses  rapidly  through  the  plasma  as  shown  by  the  B-field 
plateau.  For  z  <  z„  on  the  other  hand,  the  plasma  density  is  low,  its 
temperature  is  high  and  it  expands  in  the  axial  direction.  This  expansion 
reduces  the  magnetic  field  inside  the  plasma,  slowing  down  the  diffusion 
process  which  is  already  weaker  than  or  the  low-temperature  side  (at 
z  -  0).  The  ensuing  confinement  effect  for  z^  <  z  <  z^  is  seen  on  the 
velocity  and  density  profiles  in  Fig.  3c, d.  The  density  profile  in  Fig.  3d 
stays  nearly  a  constant  due  to  the  jxB  slowing  down  effect  in  the  above 
domain.  When  the  j  x  B  force  was  artificially  suppressed  in  the  calculations, 
no  density  plateau  was  observed. 
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As  time  increases,  j  x  B  may  be  large  enough  so  that  the  velocity  profile 
becween  z (higher  velocity)  and  2 v  (lower  velocity)  changes  significantly. 

In  that  case,  the  distance  ( z -  z^)  will  shrink  to  0  during  the  ion  beam 
pulse  duration.  The  front  velocity  Vp  will  then  exhibit  a  sudden  change  in 
value  from  v(2v)  to  vCz^)  as  (zv~  z^)  goes  to  zero.  The  whole  anode  plasma 
at  z  <  is  accelerated  during  that  time,  reaching  velocities  of  tens  of 

cm/ysec.  This  situation  occurs  especially  at  small  radii  where  both  j  and  B 
are  larger.  For  the  case  of  Fig.  3,  the  velocity  Vp(z^)  remains  approximately 
constant  at  16  cm/ysec  from  10  to  30  nsec  whereas  vp(zv)  gradually  decreases 
from  12  cm/ysec  to  5  cm/ysec.  Later  in  the  pulse,  as  (zv  -  z^)  approaches 
zero,  the  plasma  moves  at  velocities  up  to  30  cm/ysec  at  t  »  50  nsec.  At 
large  radii  however,  j  x  B  is  smaller  and  does  not  generate  a  velocity 
differential  between  2^  and  zy  large  enough  to  cause  complete  shrinkage  of 
this  gap.  The  plasma  front  velocity  Vp  is  nearly  equal  to  10  cm/ysec  and  does 
not  undergo  any  abrupt  changes.  The  plasma  front  location  as  a  function  of 
radius  is  shown  in  Fig.  4  at  various  times.  We  see  that  the  experimental 
observations  of  large  velocities  late  in  the  pulse  as  well  as  larger 
acceleration  near  the  anode  center  are  all  in  qualitative  agreement  with  the 
results  of  the  present  model.  The  results  of  the  present  model  have  been 
checked  for  their  sensitivity  to  the  following  variations. 

Initial  plasma  mass;  the  mass  was  changed  from  3.5  x  10  g/cm  to  3.5  x 
10  g/cm  for  r  »  1  cm  while  keeping  the  same  initial  density  gradient.  The 
thermal  expansion  in  the  early  stages  is  very  similar  and  resulted  in  the  same 
value  for  the  front  velocity.  Then,  in  both  cases,  the  current  density  peaks 
at  the  foot  of  the  density  gradient  and  gradually  erodes  it.  In  time,  this 
plasma  expands  and  as  mentioned  previously  is  slowed  down  by  the  jxB  force. 

'Jp  to  times  of  30  nsec,  there  are  no  major  differences  between  the  two  mass 
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I  mm 
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time.  Calculations  have  been  made  at  three  different  radii 


density  cases.  However,  later  in  the  pulse,  the  reduction  in  density  near  2  ■ 
0  for  the  more  massive  target  is  slower  and  once  the  plasma  starts  to  be 
accelerated  from  the  back,  it  reaches  lower  velocities  since  there  is  more  of 
it.  An  integral  of  the  momentum  equation  gives  the  average  velocity  in  the 
later  stages  when  the  total  pressure  differential  across  the  plasma  can  be 
neglected 


<  v  >  *  - =■  t 

inn  2  aCC 

200  it  m  r 
0 

where  tacc  is  the  time  over  which  acceleration  from  the  back  takes  place. 

This  equation  gives  also  a  very  rough  check  on  the  code  (because  of  the 
assumption  on  the  pressure  differential).  For  example,  for  the  case  of  Fig. 

3,  we  find  that  for  I  <*  1.5  MA,  tacc  *  25  nsec,  r  *  1  cm  and  mQ  *  3.5  10“^ 

2 

g/cm  ,  then  <  v  >  *  26  cm/ysec  whereas  the  code  gives  23  cm/ ysec  at  t  *  50 

2 

nsec.  For  mQ  »  8.5  x  10  g/cm  ,  the  observed  pressure  differential  is  not 
negligible  and  we  find  that  <  v  >  =  10.7  cm/ysec  whereas  the  code  gives 
15  cm/ysec  at  t  *  50  nsec. 

Maximum  temperature  in  the  low-density  plasma:  a  low-density  plasma 
-7  -i 

(p  <  10  g/cm  )  may  be  present  in  the  calculations  at  both  boundaries,  on  the 

cathode  side  naturally  and  on  the  anode  side  late  in  the  pulse  when  the  whole 

plasma  is  being  pushed  from  behind  by  the  (jxB)  force.  If  current  is  still 

allowed  to  flow  in  this  region,  the  temperature  can  reach  several  keV. 

15  -3 

3ecause  of  the  low  plasma  density  5  10  cm  )  and  the  high  current 

6  2  9 

density  (j  3  10  A/ cm  ,  see  Fig.  3b)  the  drift  velocity  (V^  ~  5  x  10  cm/sec) 

exceeds  the  thermal  velocity  for  the  electrons  and  the  plasma  becomes 

anomalously  resistive.  In  order  to  avoid  such  excessive  heating,  a  simple 

anomalous  resistivity  model  was  used  and  a  temperature  limit  was  also  imposed 
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in  that  region.  The  present  results  were  obtained  for  T  <  200  eV.  This  limit 
affects  only  the  results  for  the  lower  mass-density  case  at  r  *  1  cm.  For  all 
the  other  cases,  the  temperature  did  not  reach  100  eV  in  that  region  (either 
because  of  the  larger  mass  or  the  lower  current). 

Conclusion 

In  conclusion,  it  has  been  shown  for  the  first  time  that  anode  plasma 

expansion  can  be  explained  in  terms  of  a  magneto-hydrodynamic  model  and  that 

the  main  observations  of  the  time  history  of  the  expansion,  large  final 

velocities  (30  cm/ ysec)  and  preferred  acceleration  near  the  anode  center  are 

well  reproduced  by  the  model.  Improvements  to  the  present  model  should 

include  a  detailed  treatment  of  the  resistivity  in  the  low-density  regions 
1 5  —3 

(ng  <  10  cm  )  as  well  as  a  two-temperature  model.  The  important 
implication  of  this  work  for  plasma  diodes  is  that  there  is  a  limit  on  how 
small  high  power  diodes  can  be  made  and  how  short  the  pulse  length  should  be 
before  running  into  limitations  due  to  anode  plasma  expansion.  Other  diodes 
besides  the  Pinch-Reflex  diode  should  also  be  investigated  in  that  respect 
using  models  similar  to  that  one  presented  here. 
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